Chiral Metal as a Heisenberg Ferromagnet 



On 
On 

> 
O 

in 



I 

-a 
c 

o 
o 



> 

m 
O 



ON 
-I— > 



i 

C 

O 

o 



H. Mathur 1 ' 2 

1 Department of Physics, Case Western Reserve University, Cleveland, Ohio 44106-7079 
2 Institute for Theoretical Physics, University of California at Santa Barbara, Santa Barbara, CA 93106-4030 

(February 1, 2008) 

The two-dimensional surface of an integer quantum hall multilayer is mapped onto a Heisenberg 
spin-chain with ferromagnetic coupling. Using this mapping it is shown non-perturbatively that the 
surface states constitute a very anisotropic metal in the infinite size limit. For multilayers of finite 
size, two diffusive mesoscopic regimes are identified and the conductance fluctuations are calculated 
perturbatively for both. The Heisenberg spin-chain representation is used to study the directed 
wave problem and the exact result is obtained that the mean-square deflection of a directed wave 
grows as the square root of the propagation distance. 
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Dirty electronic systems exhibit a variety of phases 
many of which are not well-understood 0] , • For exam- 
ple electrons moving in two dimensions under the influ- 
ence of a magnetic field exhibit a sequence of localization- 
delocalization transitions as the magnetic field or some 
other system parameter is varied. Understanding of these 
transitions, which underlie the quantum hall effect, is 
based largely on numerical simulation Q. 

The purpose of this Letter is to study the two dimen- 
sional electronic states that live on the surface of an inte- 
ger quantum hall multilayer (sometimes called the bulk 
quantum hall effect). At the edge of each quantum hall 
layer the electrons circulate in one sense only and may 
be modelled as non-interacting chiral fermions j|] . If the 
layers are coupled by tunneling, the surface states com- 
prise a two-dimensional chiral electronic system (see fig 
1). The Bechgaard salts are a natural realization || and 
it is also possible to fabricate an appropriate semicon- 
ductor heterostructure || . It will be seen below that the 
chiral model used to describe the surface states of the 
quantum hall multilayer (eq 1 below) bears some simi- 
larity to a model for the quantum hall effect introduced 
by Chalker and Coddington J^J. However there is a cru- 
cial difference which makes the chiral model tractable 
whereas the Chalker-Coddington model has defied solu- 
tion. 

The key question from the point of view of quantum 
transport is whether the surface electronic wave functions 
are localized or extended along the direction of the field 
(the z— direction in fig 1). This determines whether the 
system is metallic or insulating since transport along the 
z— direction would be dominated by the surface. In cither 
case, to fully characterize the transport it is not sufficient 
to study the disorder-averaged conductance: Localized 
electronic systems generally possess a very broad distri- 
bution of conductances. Although in contrast metallic 
systems of ordinary (rather than chiral) electrons do not 
have a broad conductance distribution, the conductance 



fluctuations of finite-sized or mesoscopic metallic grains 
have remarkable universal properties || (for example, the 
typical fluctuations are of order e 2 /h— independent of the 
mean conductance or other details of the sample) . 

Quantum transport in the chiral model has previ- 
ously been studied numerically || and by mapping onto 
field theories |To|-|l^] . In particular Balents et al. fj^l 
have mapped the system onto a one-dimensional su- 
persymmetric ferromagnetic spin-chain using the well- 
established supersymmetric technique for performing dis- 
order averages . Using this mapping they are able to 
establish the important result that the two-dimensional 
chiral model is metallic (in the sense that the conduc- 
tance scales ohmically) even for arbitrarily strong dis- 
order. This should be contrasted with conventional elec- 
tronic systems in which metallic scaling is associated with 
large conductance and which are generically localized in 
two dimensions [Q. The surface states of the bulk quan- 
tum hall system are thus revealed to be a novel metallic 
phase with interesting localization and mesoscopic fluc- 
tuation physics that awaits exploration. 

The chiral model (eq 1 below) that describes the sur- 
face of a quantum hall multilayer is of interest from 
another quite distinct point of view. It should apply 
whenever waves propagate in a medium that is suffi- 
ciently anisotropic to warrant neglect of backscattering 
in one direction. The problem of waves propagating in 
an anisotropic medium has been the focus of much at- 
tention and is known as the directed wave problem (see 
for example ref and references therein) . Although 

directed waves are described by the chiral model the ques- 
tion one asks in this context is very different: The elec- 
tron is assumed to be localized at a single point initially. 
Thereafter it moves ballistically in the chiral x— direction 
and the wave function spreads (presumably diffusively) 
in the transverse z— direction. The interesting questions 
here concern the growth of the wave packet width and 
the fluctuations in the position of the centre of the wave- 
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packet (denoted [< n > 2 ]i mp ; the notation is explained 
below). The broadening of the width can be easily cal- 
culated and the answer has been known since the 1970s 
p6| . [< n > 2 ]im P is more difficult to calculate as the av- 
erage of four green's functions is now required; it has not 
been evaluated within this model previously. 

In this paper a new approach to disorder averaging is 
introduced which is distinct from the conventional replica 
or supersymmetry methods and is especially adapted to 
this system. Using this method it is possible to map the 
chiral model (in the limit of infinite size) onto a much sim- 
pler soluble one-dimensional model: an ordinary Heisen- 
berg ferromagnet. The absence of localization for the 
chiral system established by ref |P, ^0| , ^2| then emerges 
non-perturbatively as a consequence of the well-known 
quadratic dispersion of ferromagnetic magnons. The ad- 
vantages of this mapping onto an ordinary ferromagnet 
become evident when a more difficult calculation is at- 
tempted. For example [< n > 2 ]i mp can be expressed in 
terms of the matrix elements of two-magnon states of the 
ferromagnetic representation. Although magnons inter- 
act, it is not difficult to obtain the two-magnon eigen- 
states of a ferromagnet [[l7|. Carrying out such a calcu- 
lation leads to the exact result that [< n > 2 ]i m p grows as 
x 1 / 2 . This exact result agrees with the numerical simu- 
lations of [Gl| but contradicts those of ; it also agrees 
with results obtained from a simplified lattice model of 
directed wave propagation introduced by Saul, Kardar 
and Read [ p^| , pO|pl] ] . In their model special assump- 
tions are made about the disorder which make it possible 
to directly evolve the probability function (the modulus 
square of the wave function) without reference to the 
wave function itself. The results of this paper show that 
the same behaviour results for a more generic disorder 
distribution. 

Finally the properties of mesoscopic quantum hall mul- 
tilayers are studied. Due to the anisotropy of the surface 
states two diffusive mesoscopic regimes can be identi- 
fied. Following the terminology of ref p2] in the one- 
dimensional diffusive regime electrons are typically able 
to wind around the sample in the chiral direction many 
times before diffusing across in the 2— direction. The op- 
posite limit of sufficiently large circumference that wind- 
ing paths are rare is called the zero-dimensional regime. 
The mesoscopic regimes are difficult to analyse due to 
the winding paths. Spectral correlations were computed 
in ref ]l2|| in the zero-dimensional limit; however spec- 
tral correlations are difficult to probe experimentally. 
Here conductance fluctuations are calculated in the two 
mesoscopic regimes using the standard methods of im- 
purity averaged perturbation theory Universal con- 
ductance fluctuations of order e 2 /h are found in the one- 
dimensional limit coinciding with the result for ordinary 
metals. Interestingly, there is a crossover to much larger 
fluctuations in the zero-dimensional diffusive limit. 

The precise model used for the multilayer surface is 



now described. Each quantum hall layer has as many 
edge channels as there are filled Landau levels. For 
simplicity we shall focus on the case of just one filled 
Landau level. The surface of the multilayer is then a 
two dimensional chiral electronic system governed by the 
Schrodinger equation 

(-iv-^ + V n (x)-E^ G R (n,x;n',x') 

-t{G%{n + 1, x; n', x') + Gf (n - 1, x; n', x')} 

= -iv5 n>n 'S(x — x'). (1) 

Here G R is the retarded Green's function at frequency 
E. t produces interlayer hopping and V n is the disor- 
der potential. Units are chosen so that the edge velocity 
v — 1, the interlayer separation a = 1 and h = 1. The 
anisotropy and chiral character of the model are reflected 
by the fact that the equation is first order in the chiral 
x— direction whereas it is second order in the transverse 
direction. 

A mapping onto a one-dimensional problem is obtained 
by noting that in the limit of infinite size (but in that 
limit only) the circumference C — > 00 and the retarded 
Green's function obeys the chiral boundary condition 
G R (n, x; n' , x') — for x < x' . Due to this boundary 
condition, it is possible to interpret eq (1) as the time- 
dependent Schrodinger equation for a one-dimensional 
tight-binding model with the chiral co-ordinate x identi- 
fied as time and G R (n,x, ;n',x') identified as the time- 
domain retarded Green's function. Note that the on- 
site energies of the tight-binding model fluctuate in time. 
It is convenient to make a gauge transformation G R — > 
G R exp[i7„(x) - ij n > (x')] where dj n (x)/dx = V n (x) - E. 
In this gauge eq (1) becomes 

Q 

— i— G§(n, x; n', x') — t n (x)G^0n + 1, x; n', x') 
ox 

-tl_ 1 {x)G*{n-l,x\n , ,x , ) = -i5 n>n ,8{x-x'). (2) 

Heret„(a;) = texpi[j n+ i(x)— 7„(x)]. The random poten- 
tial is thus eliminated but the hopping matrix elements 
are now modulated by random phases. 

The tight-binding model can be rewritten in second- 
quantized language by introducing c R ^ and c R which cre- 
ate and annihilate fermions on site n of the tight-binding 
lattice and which evolve in the chiral time direction ac- 
cording to the Hamiltonian 

fcfdO*) = E{*«(*) C n tc n+l + (3) 

In this language the Green's function is given by 
G R (n,x;n',x') =< 0\c R Pexp (i J dxxh^x^ c£?|0 > 

(4) 
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for x > x' . The symbol Pexp() denotes a time-ordered 
exponential. Note that there is no vacuum amplitude 
in the denominator because we are calculating the single 
particle green's function rather than the propagator for 
a single particle added to a filled fermi sea, which is the 
object usually studied in many-body physics \22\. The 



absence of a denominator in eq (4) is a crucial simplifi- 
cation that allows disorder averaging as discussed below. 

The complex conjugate of the Green's function can be 
calculated from an expression analogous to eq (4) by in- 
troducing conjugate fermions, c Ajf and c A which evolve 
according to a conjugate hamiltonian (obtained from eq 
(3) by making the replacements c R — > c A , t — > — t* and 
t* -* -t). 

Localization, or its absence, is established by cal- 
culation of the disorder-averaged diffuson propagator 
\G R (n, x; n', x')\ 2 . For this purpose it is neccessary to 
simultaneously introduce both sets of fermions evolving 
according to the total hamiltonian 



h ld (x) = hf d + hf d = Y\t n A n + t* n Al}. 



(5) 



Here A n = c R ^ c R +1 — Cn+\ c n- The diffuson is then given 
by 

\G(n,x;n',x')\ 2 = 
< 0\c A c R Pexp (i f dxthwixx)) c$$\0 > . (6) 



Eq (6) provides an exact formal expression for calcu- 
lating the diffuson for a given disorder. The task now 
is to average over different realizations of the disor- 
der potential which is taken to be gaussian white noise 
with correlations [t n (x)t m (x')]i mp = DS n . m S(x — x') and 

[*n(aO]imp = [tn(x)t m (x')]i mp = 0-0, where [...Jimp 

denotes an average over disorder. Due to the absence of 
a denominator in eq (6) this average is easily performed. 
It is neccessary only to evaluate 



Pexp \ i J dxihid(xi) 



= exp{-h in t(x - x 1 )} 



where 
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A n Al) 



(8) 



Eq (8) can be verified by expanding the exponentials in 
eq (7). The rough content of eq (7-8) is that for cal- 
culating averages, the fermions may be taken to evolve 
according to an effective Hamiltonian /i mt which is not 
random and does not depend on x. It is an interacting 
hamiltonian since A n is bilinear. 

Consider a state in which a single site n is simultane- 
ously occupied by R and A fermions— below this state will 



be identified as a magnon localized at site n. The effect 
of hi n t on such a state is to cause both fermions to hop 
together onto a neighbouring site. This physics can be 
brought out clearly by defining = ^(c R ^c R ~ c n cA ^)i 
J+ ee J- + ijy ee c^c A \ J- ee (J+)t, which satisfy the 
su(2) algebra, and N n = c^c^ + c A c A ^ , which commutes 
with all the J's. In terms of these operators 
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-N n N, 



(9) 



—evidently a Heisenberg ferromagnet. As usual the vac- 
uum |0 > is the ground state and exact low- lying excita- 
tions are magnons obtained by constructing plane- waves 
from the localized magnons mentioned above. Explicitly, 
a magnon of wave- vector k is given by ^ n exp ikn\0 > 
and has eigenvalue 2D(1 — cosfc). 

The exact diffuson propagator can now be straightfor- 
wardly calculated by substituting eq (7) in eq (6) and ex- 
panding c R }c A }\0 > and < Q\c A c R in terms of magnons. 
The result is 

[|G(n,a ; ;n'^')| 2 ]i mp = ^-^)x 

/+*■ dfc 
— exp ik(n - n') exp{-2P(l - cos k){x - x')} (10) 

in agreement with eq (98) of Balents et al. | fi2[ . The phys- 
ical content of eq (10) is that the electrons move ballisti- 
cally in the chiral direction and diffuse in the transverse 
direction (diffusion constant = D); this is revealed, e.g., 
by using eq (10) to calculate the density response func- 
tion m. 

Next consider the directed wave problem. The elec- 
tron is assumed to be initially localized at the origin of 
the co-ordinate system n' = 0,x' — 0. After it moves 
ballistically in the chiral direction to a location x, the 
amplitude to be in layer n is given by G R (n, x;n' = 
0,x' — 0) and the position of the wave-packet centre 
< n >= J2 n n\G (n, x; 0, 0)| 2 . By symmetry, evidently 
[< n >]i mp = 0; and the mean-square deflection of the 
wave-packet centre is therefore 



(7) [< n > 2 ] imp = Y, nm l\G{n,x;0,0)\ 2 \G{m,x;0,0)\ 2 



(11) 



The large x asymptotic behaviour of [< n > 2 ]i mp is de- 
sired. 

To perform this calculation it is neccessary to intro- 
duce two sets each of R and A fermions (c R , c A , d R , d A ) 
which evolve according to the hamiltonian of eq (5) but 
with A n — * c^c^ +1 — c A \\C A + (c — > d). Repeating the 
previous arguments, 



[|G(n,x;0,0)| 2 |G(m,a;;0,0)| 2 



limp 



< 0\c A c R did R exp{-h int x}d R U A ^ R ^c A ^\0 



> 



(12) 
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Here fti„ t is given by eq (8) but with A n redefined as 
above. The important states are (localized) two-magnon 
states which are of two kinds. In the first type a c R — c A 
pair occupies one site while a d R — d A pair occupies an- 
other; in the second type c R is paired with d A and c A with 
d R p^ j. These states are closed under the action of h int 
which generally causes paired fcrmions to hop together to 
a neighbouring site. An exception is when the two pairs 
occupy adjacent sites, in which case the fermions may 
switch partners and a state of one kind is transformed 
into the other. Thus the two-magnon states define a sort 
of two-body problem with a contact interaction of a kind 
familiar from ordinary ferromagnetism [fL7| . Following 
the standard method explicit forms for the two-magnon 
eigenstates are obtained which can be used to straightfor- 
wardly evaluate eq (11-12); details will be given elsewhere 
pl[ . Contrary to the naive assumption that the centre 
undergoes a random walk, and hence [< n > 2 ]i mp ~ x, 
exact calculation reveals that [< n > 2 ]imp = \J {2Dx)/ir. 

Finally consider multilayers that are of finite size in 
both chiral and transverse direction thus permitting elec- 
trons to wind around the sample in the chiral direc- 
tion. During each circumnavigation, the electrons will 
typically diffuse a distance \J DC in the transverse di- 
rection. Consequently there are two distinct regimes 
depending on whether N <C \JDC (zero-dimensional 
limit) or N ^s> V DC (one-dimensional limit). Here 
N = number of layers in the multilayer. Finite sized 
samples are difficult to analyse due to the complex in- 
terference produced by winding paths. However, the 
conductance fluctuations can be calculated using dia- 
grammatic perturbation theory within the approxima- 
tion normally used for diffusive electrons The re- 
sult is [(5g) 2 ]imp — A(e 2 /h) 2 (one-dimensional limit) and 
[{Sg)%m P = A'{(CD)/N 2 }(e 2 /h) 2 (zero-dimensional 
limit) . Here A and A' are constants of order unity which 
we have not computed explicitly p5[ . Conductance fluc- 
tuations are of order e 2 /h in the one-dimensional limit 
as they are for ordinary metallic grains. The larger 
non-universal result for the zero-dimensional limit may 
be interpreted as follows: The electron moves a dis- 
tance N 2 /D <C C in the chiral direction before it dif- 
fuses into the phase-randomizing probes. The sample 
therefore breaks up into CD/N 2 incoherent blocks, each 
with independent conductance fluctuations of order e /ft. 
Although the perturbative results given here have an 
appealing physical interpretation, the validity of these 
results deserves further study via simulation or non- 
perturbative analysis (which may be possible in the zero- 
dimensional limit |2(| ) . By the familiar "ergodic hypoth- 
esis" of mesoscopic physics M we expect the statistical 
fluctuations of the conductance calculated here would be 
experimentally manifested as fluctuations in the conduc- 
tance of a given specimen when a tunable parameter, e.g., 
the magnetic field, is varied. 



It is interesting to compare the method of this pa- 
per with other field theory representations of the chi- 
ral model. Because single particle properties are be- 
ing calculated, eq (3—8) would remain valid even if the 
fermions were replaced by bosons {c R 1 c A — > b R ,b A ) 
yielding an interacting boson representation of the chiral 
model rather than the Heisenbcrg ferromagnet analysed 
here. If fcrmions and a redundant set of bosons are in- 
troduced, the supersymmetric representation of Balents 
et al. |l^] results. It is very convenient to be rid of the 
unneeded bosons, e.g., while analysing the directed wave 
problem; but it is important to emphasize that the bosons 
are optional only for the infinite system. For a finite cir- 
cumference the boundary conditions on eq (1) change and 
it becomes neccessary to introduce bosons whether oper- 
ator methods (this paper) or functional methods (Balents 
et al.) are used. 

In summary, the main results of this paper are: a map- 
ping of chiral waves onto a Heisenberg spin-chain with 
ferromagnetic coupling; calculation of the diffuson prop- 
agator, eq (10), which shows the surface of a quantum 
hall multi-layer to be an anisotropic metal in the infi- 
nite size limit; diagrammatic calculation of conductance 
fluctuations in two diffusive mesoscopic regimes for finite 
sized multilayers; and the exact result that the mean- 
square deflection of a directed wave grows as the square 
root of the distance it has propagated. 
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